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Abstract 

We analyze a supersymmetric system with four flat directions. We observe several interesting properties, 
such as the coexistence of the discrete and continuous spectrum in the same range of energies. We also 
solve numerically the classical counterpart of this system. A similar analysis is then done for an alike, 
but non-supersymmetric system. The comparison of theses classical and quantum results may serve as a 
suggestion about classical manifestations of supersymmetry. 
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Great interest has been recently put on the 
Supersymmetric Yang-Mills Quantum Mechanics 
(SYMQM), which results from a dimensional reduc- 
tion of the supersymmetric Yang-Mills field theories 
to a single point in space. The use of simple numerical 
methods gives a good understanding of the N = 2, 
D = 2 and D = 4 models Q] 12 13] ■ The ultimate 
goal of such analysis is the D = 10, SU(N — ► oo) 
model, which is conjectured to be in relation with the 
M-theory g]. 

The system addressed here was used by de Wit, 
Liischer and Nicolai in [5; 6_ as a simplified model of the 
supermembrane hamiltonian in an extensive discussion 
on stability of supermembranes, structures present in 
the M-theory. It is one of the simplest supersymmet- 
ric models with flat directions. It thus contains several 
interesting properties of SYMQM, such as the coin- 
cidence of the discrete and continuous spectrum. The 
aim of this paper is to analyze in detail this system with 
the method proposed by Wosiek pQ and to compare 
it with an analogous, but non-supersymmetric model. 
The latter was also already described in the literature 
[7], because of its interesting feature: all its eigenstates 
are localized, even though the potential is zero on an 
unbounded set. 

The paper is composed as follows. We start by intro- 
ducing the quantum systems: the non-supersymmetric 
one and the supersymmetric one, and we continue by 
presenting their classical counterparts. The compari- 
son in the quantum regime is based on the analysis of 
numerical spectra, therefore we proceed by describing 
the method used for calculating spectra, and then give 
a detailed study of all symmetries present in theses 
systems in order to fully understand the degeneracies 
which may appear. Consequently we discuss the re- 
sults of this comparison. Finally, we turn to the classi- 
cal regime and analyze the classical bosonic and super- 
symmetric trajectories. Comparing them enables us to 
rediscover the differences between quantum systems on 
the classical level. 

Description of the systems 

In this chapter we introduce the hamiltonians of the 
considered models; first the quantum, then the classical 
one. In the following, the non-supersymmetric system 



will be called bosonic, because it contains only bosonic 
degrees of freedom. 

Quantum systems 

The quantum bosonic hamiltonian has the form 

H bo sonic=pl+Pl+X 2 y 2 . (1) 

The potential x 2 y 2 has four flat directions. Usually, 
we expect that systems with potentials equal to zero 
on an unbounded set, have continuous spectrum. Nev- 
ertheless, in our case, the spectrum turns out to be 
exclusively discrete due to the quantum fluctuations 
in the transverse directions. Let's consider a particle 
moving in one of the valleys, say x > 0. The trans- 
verse potential is the potential of a harmonic oscillator 
with an equilibrium position at y = and frequency 
proportional to the distance from the center ui ~ \x\. 
The zero-mode energy of such fluctuations increases 
linearly with |a;|, when the particle is moving deep into 
the valley. Therefore, the particle is exposed to an ef- 
fective potential barrier, which prevents it from escap- 
ing. We will now recall [7J the solution of the station- 
ary Schrodinger equation in the Born-Oppenheimer ap- 
proximation in order to explicitly show the above. The 
energy of a quantum harmonic oscillator of frequency 
| a; | is equal 

E x = \x\(n y + -). 

The wavefunction can be written as a product of two 
functions 

$(x,y) = a(x)f3 Xtny (y), 

where a(x) accounts for the onward motion, and 
Px,n y (y) describes the transverse harmonic fluctuations 

d 2 1 

(-Qy2 + X V )Px,n y (y) = 2\ X \{ny + -)P x ,n v (y)- 

Thus, (3 x ,n y (y) is a well-know eigenfunction of quantum 
harmonic oscillator 

P x ,n y (v) = (V^2 n «n y l)-^ 2 H ny (y/2\x]y) e"^'. 
a(x) satisfies a stationary Schrodinger equation 

<9 2 1 

(~g^2 + 2 \ X \( n V + l}))^) = E n^n y a(x). 
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We now proceed independently for x > i x < 0. Let 
us introduce a new variable z 

* = <- + 5> I/3 (*-2^l)' 

The equation for a(x) can be now rewritten in this 
variable as 

a{z) — 2za{z) = 0. 

This is the Airy equation. Its solutions are even and 
odd Airy's functions. In order to obtain a full solution 
we have to match together solutions from x > and 
x < valleys. This gives us the condition of energy 
quantization, which turns out to be the equation for 
the zeros of Airy's function 

This condition has a straightforward interpretation. 
The spectrum is discrete, because the Airy's functions 
are analytic and therefore have countable many zeros 
on the real axis. The few first approximate eigenen- 
ergies are shown in table [5] on page [H together with 
some exact ones. The compatibility is within 7%. The 
assumption, that the particle moves only in one of the 
valleys, breaks the space symmetries of the system, so 
the approximated energies will not have degeneracies. 

The quantum supersymmetric hamiltonian has the 
form 

H susy =p 2 x +p 2 y + x 2 f + (x + iy)f + (x - iy)f, (2) 

where and / are the fermionic creation and annihila- 
tion operators respectively. The operator f'f does not 
commute with the hamiltonian ^ , so the fermionic oc- 
cupation number is not a good quantum number and 
cannot be used to label the eigenstates. This gener- 
ates non- vanishing, non-diagonal matrix elements in 
the fermionic occupation number representation. 
There exists one supersymmetry generator 

Q = Q f - -4((p* + ipy)f + (p x - ip y )P + xy[f\ /]), 



such that 

H = {Q\Q}, 

and 

[H,Q} = 0. 

Although the supersymmetry concerns particles - 
bosons and fermions, it is sometimes useful to think 
about the fermionic degree of freedom as an equivalent 
spin projection on the OZ axis. In this language the 
hamiltonian and the supersymmetry generator take the 
form 

Hsusy =Pl+Pl + x 2 y 2 + a x x - <r y y, 

Q = & = -^=(p x a x +p y (Jy +£ya z ). 

The analytic analysis of the system ([2]) is difficult be- 
cause of the mentioned non-diagonal matrix elements 
in the fermionic occupation number representation. 
We can nevertheless extract some useful information 
by analyzing a simpler and more regular model. It 
will enable us to demonstrate the effect of coexistence 
of the discrete and continuous spectrum in the same 
range of energies. Let's consider the hamiltonian 

H model =p 2 x +p 2 y + x 2 y 2 + x[f^f]. (4) 

We will use now, as in the purely bosonic case, the 
Born-Oppenheimer approximation. It has to be noted 
that although the model hamiltonian (j4]) does not have 
the space symmetries present in the original system (|2|). 
it does not spoil its usefulness, because the approxima- 
tion breaks these symmetries anyway. Let's consider 
the motion in one of the valleys, say x > 0. The trans- 
verse potential is a potential of a quantum supersym- 
metric harmonic oscillator of frequency u> ~ |x|. The 
hamiltonian of such oscillator is 8J 

1 2 1 2 2 |l| j 

H SUS y oscillator ~Q^y 9 ^ ' 

Its energies depend on one quantum number n and are 
equal 

E n = \x\n. 

Similarly to the bosonic case, we write the wavefunc- 
tion as a product of two functions 

*(x,y) = a(x)(3 x . n (y), 
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where a(x) accounts for the onward motion and (3 x ,n{y) 
fulfills the equation 

(-^2 + x 2 y 2 - \x\[f, /*])&,»(») = 2n\x\[3 x , n (y). 

This gives us the equation on a(x) 
d 2 

\~~q~2 + 2n\x\)a n {x) = Ea n (x). (5) 

From the form of ([5]) we can conclude that the spec- 
trum has two coexisting parts: a continuous part for 
n = and a discrete one for n > 1. The zero energy 
of the zero-modes enables the particle to penetrate the 
valley to any depth. The same effect is present in the 
original system ([2|), so it seems that ((4J) is a good can- 
didate to model this feature. 

A detailed study of spectra of both bosonic and su- 
persymmetric quantum models is described later on. 

Classical systems 

The classical bosonic hamiltonian has the form 

H boz =pl +p 2 y + x 2 y 2 . (6) 
We easily find the equations of motion 

x = 2p x , 

V = 2Py, (?) 
p' x = -2xy 2 , 

Py = -2yx 2 . 

These equations imply that particles which move ex- 
actly along the coordinate axes have a constant mo- 
mentum component along these axes. This, of course, 
corresponds to a free motion in these directions. We 
will show later that any other trajectory is bounded 
in the sense that it always returns to the center of the 
potential. This property of the x 2 y 2 potential can be 
easily investigated using the so-called hyperbolic bil- 
liard. The authors of the article [H] have considered a 
classical, free, spinless particle moving as a free particle 
in a bounded set 

D = {(x,y)\x >0Ay>0Ay< l/x}, 



which is an idealized version of the potential in 
They showed that such particle cannot escape through 
the valleys, it always turns back. They also proved 
that it is possible to systematically find all closed or- 
bits. Every orbit can be unambiguously labeled by a 
set of symbols which correspond to reflection from par- 
ticular walls of the potential. Such a list of all orbits is 
particularly useful in the semiclassical calculations by 
Feynman integrals. 

The classical equations of motion of the supersym- 
metric system can be obtained from the quantum 
hamiltonian @ using the Ehrenfest's equation for 
mean values 

x = 2p Xl 

y = 2p v , 

Px = -2xy 2 ~2S x , 

p y = -2yx 2 + 2S yi (8) 

S x = -2yS z , 

Sy = 2xS Z7 

S z = 2xS y + 2yS x . 

Equations ((5|) can also be derived in a different way. 
As was shown in the paper by Berezin [10] . for the 
description of classical dynamics of a nonrelativistic 
particle with spin ^ one needs an enriched phase space. 
This description uses three, anticommuting, dynamical 
variables which belong to a Grassmann algebra G3 
with three generators. 

&&+£i& = 0, k,l= 1,2,3 

In particular £ 2 = 0. We can use relations between 
quantum spin projection operators with their classical 
counterparts expressed in terms of £. We thus obtain 
a classical form of the supersymmetric hamiltonian 

H SUS y =P 2 c +pI + x 2 y 2 - i(xe xlm - ye y i m )^ m - 

The time dependence of the dynamical variables 
can be derived from the Hamilton equations : 

§ = f> # = -ff> t = -^re q = x,y 

and p = p x ,p y . The appearance of the imaginary i 
in the classical equations is inherently related to the 
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nature of the Grassmann variables. 



X 


= 2p x , 




y 


= 2 P y, 




Pi- 


= -2xy 2 




Py 


= -2ya; 2 




ii 


= -2y&, 






= -2^ 3 , 
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= 2x& + 


2^i- 



(9) 



It is useful to replace the Grassmann variables in 
these equations by the spin projections Sk, which are 
the physical observables. 
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Figure 1: Vector field V about which the spin rotates. 



We get rid of the time derivatives using the equa- 
tions of motion © . The equations for Sk close and we 
get the already obtained set of equations ([5]). In this 
language, the hamiltonian takes the form 



H suay — p x 



2 , 2 2 

■ P v + x y 



(10) 



Furthermore, equation (|10p can be rewritten in a form 
where the spin precession is evident. To this end we 
define a vector field, which is space-dependent 



Then 



V=(V x ,V y ,V z ) = (2x,-2y,0). 



Hsusy =pl+P 2 y + x 2 y 2 + S-V, 



which manifestly shows the precession effect. It is note- 
worthy that such a field is tangent to the contour line 
of the potential x 2 y 2 , see figure [TJ 

Quantum mechanics on PC 

A numerical analysis of quantum systems can be eas- 
ily carried out in an eigenbasis of occupation number 
operators, a so-called Fock basis [TJ. The occupation 
number operator can be written as (va, where a* i a 



are standard bosonic creation and annihilation opera- 
tors respectively. They fulfill well-known commutation 
relations 



[dp, a q ] = [at a\] = 0, [op, oj] = 5 p 



(11) 



p and q in (|11| are indices which label the 
bosonic degrees of freedom. Particularly, the non- 
supersymmetric system ([T]) has two bosonic degrees of 
freedom, labeled x i y. Thus, a basis state is described 
by two integers 



to, n >= 



1 



I ml 



(4r(4)"|0>, m,n>0. (12) 



Of course, in a numerical analysis it is impossible to use 
an infinite basis, so a cut-off, N cut , is needed, which 
limits the number of basis states. There are many 
ways to introduce such a cut-off. One of them, called 
a square cut-off, is to limit independently the maxi- 
mal number of occupation for each of two degrees of 
freedom by ^N cut . 

The momentum and position operators can be ex- 
pressed by creation and annihilation operators in the 
following way 



1 



V2 



(a x + a x ) 



y = -/=K + 4) 



Px 

Pv 



7^=(a x -<4), (13) 
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The action of the hamiltonian, which is an opera- 
tor function of (fT3|) , is straightforward in the Fock ba- 
sis. We can easily calculate its matrix elements. The 
eigenenergies of a quantum system belong to a set of 
eigenvalues of the hamiltonian matrix, and the eigen- 
states are the eigenvectors of this matrix. 

In the supersymmctric case we must introduce 
fcrmionic creation and annihilation operators, / and 
/t, which fulfill the anticommutation relations 



{f P J,} = {flfh = o, {U,fi} = s s 



I'll' 



(14) 



where p and q are indices which describe the fermionic 
degrees of freedom. In general, in order to ensure the 
relations (fT4")) one uses a construction by Jordan and 
Wigner [11]. In a case of only one fermionic degree 
of freedom it is not necessary. The fermionic occupa- 
tion number operator f'f, with two eigenvalues Oil, 
permit to label the basis states by a third quantum 
number 



|m, n, k >- 



\fn\m\ 



(oi) m (4) n (/ t )*io >, 



m, n > 0, k = 0, 1. 



(15) 



The above method allows also to verify the relia- 
bility of numerical results. The rate of convergence 
of eigenenergies with an increasing N cu t is a simple 
criterion. For example, we can assume that the con- 
vergence is reached when a relative change of the en- 
ergy between consecutive N cut is less than 1% of its 
absolute value. For the bosonic ground state this hap- 
pens for N cut > 100. The excited states have a worse 
convergence so all the following results were obtained 
for N cu t = 400. The dependence of the energy of the 
bosonic ground state on the cut-off is shown in figure 
[2] with a solid line. 

This dependence provides us also with an additional 
information on the character of the quantum state. 
Particularly, as it will be described later, it enables 
one to decide whether the state belongs to the discrete 
spectrum or to the continuous one. Another criterion 
of reliability is supplied by the spatial probability dis- 
tributions. One can assume that the convergence is 
reached when for consecutive cut-offs, N cu t and N' cut , 
the probability P{x) of finding the particle at the point 
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Figure 2: The dependence of the ground state energy 
on the cut-off for the bosonic system. 
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Figure 3: The dependence of the space distribution of 
the bosonic ground state on the cut-off N cu t, for N cut 
equal 16 (highest curve), 36, 64, 128 (lowest curve). 



x fulfill the inequality 

\PN aut (x) 



max 

x e i 



PN cut (x) 



< e, 



where / is an interval on which the distributions are 
calculated. Figure [3] shows a cross-section along the 
OX axis of such distribution for the bosonic ground 



6 



state for different cut-offs. The convergence is reached 
for N cut = 128. The free parameters were chosen as 
/ = [—5,5] and e = 0.01. Obviously these values are 
only exemplary and can be tuned at will. 

Symmetries of the systems 

Bosonic case and the C 4v group 

The hamiltonian of the quantum bosonic system JTJ) is 
invariant with respect to a point group of symmetries, 
which consists of a 4-fold principal axis of symmetry 
and four reflection axes perpendicular to the principal 
axis, and is called the C 4v group. It has eight elements 
e - identity, 
tt x - reflection with respect 

to the x axis , 
TTy - reflection with respect 
to the y axis , 
7T x —y - reflection with respect 
to the x = y line , 
ir x= - y - reflection with respect 
to the x = —y line , 
Ril = R - rotation through 

the angle of +-| radians, 
R n = R 2 - rotation through 

the angle of +tt radians, 
i?3^ = R 3 - rotation through 

the angle of +^ radians. 

We call a conjugacy class a subset of all elements of 
a group G that commute with all elements of G. The 
group C^ v has five conjugacy classes 



id 
2C 4 
C 2 

2S d 



{e}, 

{R, R 3 }, 

{R 2 l 

{nx,TTy}, 

\^x— yi ^ x — — y\ • 



It is known [12] , that the number of irreducibles repre- 
sentations of a group is equal to the number of classes 
of this group. At the same time, the sum of squares 
of dimensions of this representations must be equal to 
the order of the group. All this implies that C 4v must 
have four 1-dimensional and one 2-dimensional irre- 
ducible representations. The defining representation is 



2-dimensional and is explicitly shown in p^|) . Each 
group has also a 1-dimensional, trivial representation 
composed of ones. 



1 
1 



R 2 



-1 






-1 



R- 
R 3 



-1 
1 

1 

1 



7T X — — y — 



1 

-1 

-1 

1 

1 

-1 

-1 
-1 



(16) 



In order to find the remaining three 1-dimensional rep- 
resentations one uses a couple of observations. First, 
every representation of the identity element has to be 
equal 1. Second, except the elements of the 2C4 class, 
all elements squared give the identity element. Thus, 
their representations must be equal +1 or —1. More- 
over, one can obtain the elements of the C 2 class by 
combining the elements of the classes 2C4, 2S V and 
2Sd (i.e. reflection with respect to the OX axis and a 
reflection with respect to the OY axis give a rotation 
through an angle of tt). So the representations of the 
C2 class must be equal to +1. As the representations 
of the 2C4 class are concerned, they can be equal +1 
or —1, because when squared they must give the ele- 
ments of the C2 class. In order that the multiplication 
table remains unchanged the representations equal to 
— 1 must appear in pairs. All these observations permit 
to identify all 1-dimensional representations of the Ca v 
group (table QJ. The eigenstates transform according 





e 


2C* 4 


c 2 


1S V 


2S d 


Ax 


1 


1 


1 


1 


1 


A 2 


1 


1 


1 


-1 


-1 


B x 


1 


-1 


1 


1 


-1 


B 2 


1 


-1 


1 


-1 


1 



Table 1: Four 1-dimensional irreducible representa- 
tions of the C 4v group. 



to one of the irreducible representation. We have four 
1-dimcnsional representations, Ai,A 2 ,Bi and B 2 , so 
we will have four series of nondegenerate states. States 
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transforming according to the representation E will be 
doubly degenerate. A state symmetric with respect to 
three reflections ir x , ir y and ir x=y , will transform ac- 
cording to the trivial representation A\. A state sym- 
metric with respect to ir x and ir y reflections but anti- 
symmetric with respect to ir x=y will transform accord- 
ing to the Bi representations. Table [2] contains first 
twelve states, their symmetries and the representation 
they belong to. 



Table 
tions. 

















e 




c 2 


2C A 


2S d 


A 


1 


1 


1 


1 


1 


B 


1 


-1 


1 


1 


-1 


C 


1 


1 


1 


-1 


-1 


D 


1 


-1 


1 


-f 


1 



3: Four 1-dimensional irreducible representa- 



Supersymmetric case 

The quantum supersymmetric hamiltonian remains 
unchanged under action of transformations which form 
the following group 

e - identity , 





x - 




y — 


y, 




f- 


■-P, 


/t- 


+ -f, 


TTy - 


x — 


> X, 


y - 


► -y, 




/-» 


p, 


P- 


> f, 


Hit 


x - 




y - 


> -y, 




f- 


♦ -/> 


P- 


>-P, 


7T x — y 


x — 


* y, 


y -> 


X, 




f- 


>ip, 


P- 


>-if, 


^x— — y 


x - 


* -y, 


y - 


<■ —x, 




/-> 




P- 


-►</, 


RlL - 

2 


x - 


-» -y, 


y -» 


X, 




f- 






-iP 


R 3 « - 

2 


x — 


* y, 


y - 






f- 


►-«/, 


P- 


+ ip. 



The aim of this notation is to highlight the similarity 
of the above group and the C± v group. In this spirit 
we call this group the C Av group. We again find five 
conjugacy classes 

id = {e}, 
2C 4 = {%%}, 

Cl = {i?7r}, 
2S V = {7T x ,TT y }, 

We can perform a similar analysis as in the bosonic case 
and obtain four 1-dimensional and one 2-dimensional 
irreducible representations (table [3] and formula H7|> . 




Since all elements of the C^ v group commute with the 
hamiltonian, one can choose one of these symmetries in 
order to define an additional quantum number. Let's 
choose the following 

TT y :x^x, y->-y, f-*P, P -> f. 

The action of n y on basis vectors is given by 

TT y \n x ,n y ,nf >= (-)"* \n x , n y , 1 - n/ > . 

One can define even and odd states with respect to n y 

\n x ,n y ,± >= -^=(|n a ,%,0 > ±(-)"« \n x , n y , 1 > . 

In such a basis the hamiltonian appears as a block di- 
agonal matrix with two sectors, and it is possible to 
diagonalize each block independently. Although we 
have introduced a squared cut-off the basis remains 
invariant with respect to the symmetries from the C Av 
group. The spectrum of the supersymmetric system @ 
will have degeneracies due to these symmetries for each 
cut-off. It turns out that all states transform accord- 
ing to the 2-dimensional irreducible representation, so 
they will all be doubly degenerate. Table[3]shows eight 
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u.p. 


unci 


B-0 


n x 


7T 

"V 


7T 

'* x— y 


" x — — y 


R 


R 2 


R A 


f a firm 


1 


1.1082 


1 1 737 

1-11 iJ 1 


_L 
I 


_L 
1 


_|_ 


_L_ 




_|_ 


_|_ 


A, 


2 


2.3788 


2.3381 




+ 


o 










o 


E 


3 


2.3788 


2.4414 


+ 


















E 


4 


3.0574 


3.2711 


+ 


+ 








+ 




Bi 


5 


3.5229 


3.4319 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


A X 


6 


4.1100 


4.0878 


+ 


















E 


7 


4.1100 


4.2950 




+ 
















E 


8 


4.8210 


4.8307 


+ 


+ 








+ 




Bi 


9 


5.0113 


4.8635 






+ 


+ 




+ 




B 2 


10 


5.1120 


5.0783 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


At 


11 


5.6947 


5.5206 


+ 


















E 


12 


5.6947 


5.8053 




+ 
















E 



Table 2: First twelve eigenstates of the quantum bosonic system with a corresponding representation of the 
symmetry group C^ v . In the column entitled 'Energy' the exact energies are shown, whereas the 'B-O' column 
contains values calculated in the Born-Oppenheimer approximation. In the columns from 4 to 10 are shown the 
parity or a lack of a given symmetry for each state. 





N cut = 72 


N cut = 128 


N cut = 200 


1 


0.2871 


0.2469 


0.1788 


2 


0.2871 


0.2469 


0.1788 


3 


1.2172 


0.8573 


0.7553 


4 


1.2172 


0.8573 


0.7553 


5 


2.2215 


1.8932 


1.5014 


6 


2.2215 


1.8932 


1.5014 


7 


3.1867 


2.7090 


2.4616 


8 


3.1867 


2.7090 


2.4616 



Table 4: Eight first energies for three different cut-offs 
for the supersymmetric case. The degeneracy due to 
the point symmetries is exact for each cut-off. 



first energies for three different cut-offs. 
Supersymmetry 

Supersymmetry can manifest itself in a given quan- 
tum system through the set of supersymmetry genera- 
tors, which form a well-defined algebra, or, on a more 
experimental level, through a characteristic structure 
(degeneracies) of the spectrum. In this paper we as- 



sume that the existence of a single ground state with 
zero energy and supersymmetric doublets of higher en- 
ergies guarantees the presence of supersymmetry. In- 
deed, in the case of system ([2]) the generators Q and 
Q* do not fulfill all relations of a supersymmetry alge- 
bra, namely Q 2 ^ 0. Moreover, Q does not conserve 
the bosonic occupation number, so the introduction of 
the cut-off breaks the supersymmetry and destroys its 
fingerprints on the spectrum. One expects its restora- 
tion in the limit N cut — > oo. However, there exist a 
way of establishing the supersymmetry for each finite 
N cut . One can remark, that, since Q is hermitian, the 
hamiltonian ([2|) is a product of generators Q specified 
for an infinite cut-off, and then cut to the desired di- 
mensions. On the other hand, one can conceive the 
hamiltonian matrix as a product of matrices of already 
cut generators Q. 



[H s 



[Q finite N cut ] [Q finite N c . 



(18) 



Of course, the hamiltonians ([2|) and (fT8|) agree in the 
limit of infinite cut-off. It turns out that the spectrum 
of (fTS)) is fully supersymmetric. A sample of results for 
three different cut-offs are shown in table [5] One can 
observe a double degeneracy due to the point symme- 
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tries of the system, as well as a supersymmetric degen- 
eracy. Moreover, one can notice in table [5] yet another 



N cut = 50 


N cut = 200 


N cut = 450 


0.0000 


0.0165 


0.0000 


0.0000 


0.0165 


0.0000 




0.0165 






0.0165 




0.3371 


0.3607 


0.0862 


0.3371 


0.3607 


0.0862 


0.3371 


0.3607 


0.0862 


0.3371 


0.3607 


0.0862 


1.0163 


0.8867 


0.3760 


1.0163 


0.8867 


0.3760 


1.0163 


0.8867 


0.3760 


1.0163 


0.8867 


0.3760 



Table 5: Twelve first eigenvalues of the matrix [H susy ] 
for three different cut-offs. The degeneracies have two 
origins: supersymmetry and point symmetries of the 
group Ci V . For cut-off, for which N cu t/2 is even, the 
supersymmetry is broken and the spectrum does not 
have the zero energy ground state. 

way of supersymmetry breaking. For square cut-offs 
with N cut /2 even, the spectrum can not contain a sin- 
gle ground state and N cut /4 — 1 supersymmetric dou- 
blets. In such situation the supersymmetry is broken 
and the ground state disappears. 

Quantum systems 

In this part we will compare the numerical spectra of 
the bosonic and supersymmetric quantum system. 

A part of the spectrum of the bosonic system (p} 
is shown in figured! The dependence of twelve first 
eigenenergies on the cut-off is depicted. Numerical val- 
ues are presented in table [2] on page [9] Consequently 
to the discussion of the symmetries of the system, the 
spectrum contains four nondegenerate series of states 
and one doubly degenerate. The degeneracies due to 
the point symmetries are present for each N cut , be- 
cause the cut-off does not spoil them. After the paper 
[13] the dependence of the eigenenergy on the cut-off 
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Figure 4: The dependence of the twelve first eigenen- 
ergies on the cut-off for the bosonic case. 
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Figure 5: The dependence of the bosonic ground state 
energy on the cut-off with fitted curves: an exponential 
one and an inverse square root. 



enables one to decide whether the state belongs to the 
discrete spectrum or to the continuous one. If the de- 
pendence on N cut is of the type — or slower then 
the corresponding state is a nonlocalized state from 
the continuous spectrum. Contrary, if the dependence 
is fast, for example exponential ~ e -N cut ^ t nen t ne 
state is localized. Figure [5] shows the dependence of 
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200 600 1000 1400 
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Figure 6: The dependence of the quantum virial on 
cut-off for ten first bosonic states. The convergence for 
the ground state is the fastest. The most distant from 
zero curve corresponds to the tenth eigenstate. 



the bosonic ground state on the cut-off and two fit- 
ted curves: an exponent f(N cut ) = a + b e~ c Naut and 
an inverse square root g(N cut ) = a + b/\/N cut . The 
first one fits much better which means that the ground 
state is localized. It turns out that all other states 
analyzed with this method belong also to the discrete 
spectrum. This is also confirmed by an analysis of the 
virial, which is a scalar product of vectors of momen- 
tum and position w — p ■ x. In the classical regime, 
the mean time derivative of a virial tends to zero on 
a bounded trajectory and explodes on an unbounded 
one. On the quantum level, the evaluation of the ex- 
pectation value of a corresponding operator in a given 
state can supply similar type of information. In the 
bosonic case one is left with the following expression 

< w bosomc >= 2((pl) + {pi) - 2{x 2 f)). (19) 

The results for ten first states are shown in figure [6] as 
a function of cut-off. We notice that all values con- 
verge to zero, which confirms that the bosonic states 
are localized. 

On the other hand, the supersymmetric spectrum 
is shown in figure [7] According to the symmetries of 
the system, all states are doubly degenerate due to the 



100 200 300 400 500 
Basis size 

Figure 7: The dependence of the thirty two first states 
on the cut-off for the supersymmetric case. 



point symmetries from the C^ v group. The supersym- 
metric degeneracy can be found in the infinite cut-off 
limit. The fitted cut-off dependence of the supersym- 
metric ground state is shown in figure [8] and suggest 
that it is a nonlocalized state. The analysis of higher 
states also confirms that they belong to the continu- 
ous spectrum. Thus, it seems that the spectrum of 
the quantum supersymmetric system is composed of 
purely nonlocalized states. But, the approximate so- 
lution of the supersymmetric model ((4} indicates that 
among them one might find localized states from the 
discrete spectrum. And indeed we have found such 
states in the numerical spectrum of the exact hamil- 
tonian ([2]). However, as their convergence is exponen- 
tial we may expect that they should cross states with 
slower cut-off dependence. The rules of quantum me- 
chanics claim that when a parameter in a hamiltonian 
is changed, the states of the same symmetry cannot 
cross. In our case, such a parameter is our cut-off, 
and as was mentioned, all states transform according 
to the 2-dimensional irreducible representation of the 
Ci v group. So these localized states must have an in- 
teresting realization. One can see them in figure [7] as 
deformations of the energy dependence on the cut-off. 
More precisely, there exist states which on some inter- 
val N cut have a constant energy (i.e. the deformations 
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Figure 8: The dependence of the supersymmetric 
ground state energy on the cut-off with fitted curves: 
an exponential one and an inverse square root. 



Figure 9: The dependence of the quantum virial on the 
cut-off for some selected supersymmetric states (first, 
second, fifth and sixth). 



for energies around E — 5). If the energy of such a 
state remained unchanged till some cut-off N cu t, then 
for N cut > N cu t it starts decreasing as j^—- However, 
from N cu t the energy of the higher, neighboring state 
assumes this value (here E = 5) and remains constant 
on some adjacent interval. In this way localized states 
coexist with nonlocalized ones in the same range of en- 
ergy. These results can again be confirmed by the virial 
analysis. We have to evaluate the operator 

< w susy >= 2(pl) + 2{f y ) - A(x 2 f) - (a x x) + (a y y). 

(20) 

In figure [5] the time derivative of virial for four su- 
persymmetric states is depicted as a function of the 
cut-off. For a majority of them the obtained values are 
far from 0. These states are nonlocalized. However, 
as was observed there can exist also some localized 
states. Taking into account their realization through 
many nonlocalized states, one concludes that the states 
for which the quantum virial has a value near for a 
given cut-off, are just the states that form the localized 
state. 

To summarize, the main observation is that bosonic 
and supersymmetric eigenstates have a different char- 
acter. The first one are localized and thus belong to the 
discrete spectrum, whereas the supersymmetric spec- 



trum is composed from both types of states. This dis- 
parity is essential in the following discussion. 



Classical systems 

The differences in the behavior of classical bosonic and 
supersymmetric trajectories are striking at the first 
glance. Similarly as in quantum regime the bosonic 
states were localized and the supersymmetric were non- 
localized, the bosonic trajectories are bounded and the 
supersymmetric escape through the flat directions of 
the potential. 

More precisely, the numerical solutions of bosonic 
equations of motions show that all investigated tra- 
jectories are bounded. They all come back to the center 
of the potential. Similarly, as in the case of the hyper- 
bolic billiard potential, mentioned earlier, only parti- 
cles that move along the coordinate axes can escape 
through the valleys. In the quantum regime, we used 
the virial theorem to decide whether quantum states 
were localized or not. Now, we would like to do the 
same using the classical virial theorem. It implies that 
for a bounded motion (with finite momentum and po- 
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Figure 10: An exemplary bosonic classical trajectory and the virial dependence on time. The mean time 
derivative of the virial has small values indicating that the trajectory is bounded. 



sition) in a potential that vanishes as - one has 

dw — — , . 

= — = U + 2K, (21) 

where w is the virial, 4i is the time derivative, - 
stands for the time averaging, U is the potential en- 
ergy and K is the kinetic energy. In order to verify 
the character of the classical trajectories, one needs to 
evaluate the quantity (f2~Tj) on a given trajectory. A 
representative bosonic trajectory and the correspond- 
ing plot of the virial are shown in figure 1101 The mean 
time derivative of the virial remains near zero, which 
indicates that the trajectory is bounded. Moreover, 
using the algorithm proposed by Schmelchcr and Di- 
akonos [14] we found some closed orbits, one of which 
is shown in figure [11] Investigation of eigenvalues of 
the Jacobian shows that this orbit is unstable. It is 
worth noting that similar research was made by Dal- 
hqvist and Russberg [15] who found a stable island in 
the phase space of the x 2 y 2 potential, therefore show- 
ing that the latter is not ergodic as was believed. 

On the contrary, all supersymmetric trajectories es- 
cape through the valleys of potential. However, there 
exist some of them which for a large amount of time 
remain in the center of potential. When the depen- 
dence of the number of returns of some trajectory on its 



energy is investigated, one finds sharp peaks suggest- 
ing that there exist trajectories with a huge number 
of returns. Again, we found one unstable closed orbit 
using the Schmelcher's algorithm, shown in figure 1111 
The virial theorem confirms the unbounded character 
of supersymmetric trajectories. The evaluation of ([2"T]) 
for a given trajectory shows that the value of the time 
derivative of the virial explodes as the particle escapes 
through the valley. Figure fT2l depicts theses results. 



The findings obtained in this part show that the be- 
havior of classical trajectories and quantum states is 
analogous. For the bosonic system quantum states are 
localized and classical trajectories are bounded. On 
the contrary, for the supersymmetric case, the major- 
ity of states are nonlocalized. The corresponding clas- 
sical motion is a trajectory which escapes through the 
valleys. The localized supersymmetric quantum states 
are accompanied by closed classical orbits. So, the dif- 
ferences observed between quantum systems are also 
found on the classical level. This pattern might indi- 
cate the existence of classical manifestation of super- 
symmetry. 
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Figure 11: Closed bosonic and supersymmetric orbits. 
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Figure 12: An exemplary supersymmetric classical trajectory and the virial dependence on time. The mean 
time derivative of the virial explodes when the particle escapes through the valley. 



Conclusion 

In this paper we investigated the properties of a system 
proposed by de Wit, Luscher and Nicolai, in both quan- 
tum and classical regimes. By comparing it with the 
bosonic system we observed several interesting facts. 

We started by introducing the quantum hamiltoni- 
ans of the bosonic and supersymmetric systems with 
the potential with four flat valleys. Using the Born- 



Oppenheimer approximation we succeeded in solving 
both systems. We concluded that a bosonic quantum 
particle moving in one of the valleys is exposed to an 
effective potential barrier which prevents it from escap- 
ing. Thus, all states are localized. Contrary, the su- 
persymmetric quantum particle can enter at any depth 
into the valleys. The supersymmetric spectrum con- 
sists of states from the continuous spectrum and, co- 
existing with them, some localized states. Then, we 
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introduced the hamiltonians of the classical counter- 
parts of these systems. For the supersymmetric sys- 
tem we obtained the equations of motion which de- 
scribe a particle with spin which precess around some 
vector field. The interesting feature is that this field 
is tangent to the contour lines of the potential. Next, 
we described the cut-off method and some criteria for 
testing the reliability of numerical results. The sym- 
metries of considered systems, discussed subsequently, 
fully explain all degeneracies of exact energies. The 
comparison of the bosonic and supersymmetric spectra 
confirmed the hypothesis based on the approximated 
solutions on the character of quantum states. We ob- 
served the noteworthy realization of the supersymmet- 
ric localized states as deformations of the energy de- 
pendence on the cut-off. We then turned our attention 
to the classical regime. We descibed the behavior of 
trajectories of the classical bosonic and supersymmet- 
ric systems. It turned out that the character of the 
trajectories corresponds to the character of quantum 
states. Bosonic trajectories are bounded, whereas the 
supersymmetric ones escape through the valleys of the 
potential which can be due to the spin precession. The 
above results show that the differences between quan- 
tum states are also present on the classical level. 

The supersymmetric potentials with flat directions 
appear also in some much more sophisticated physical 
problems, i.e. in systems connected with Supersym- 
meric Yang-Mills Quantum Mechanics (SYMQM). The 
investigated system is one of the simplest models pos- 
sessing many of the properties of SYMQM. A deeper 
analysis can help to understand, for example, the co- 
existence of the continuous and discrete spectrum in 
theses systems. Moreover, the surprising effect of spin 
precession on the motion of a classical, supersymmet- 
ric particle which enables it to escape through the flat 
directions deserves a further attention. 
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